Matter instability in modified gravity 
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The Dolgov-Kawasaki instability discovered in the matter sector of the modified gravity scenario 
incorporating a —p 4 /R correction to Einstein gravity is studied in general f(R) theories. A stability 
condition is found in the metric version of these theories to help ruling out models that are unviable 
from the theoretical point of view. For example, the theories f(R) = R + ap 2 ^ n+1 ' /R n , where a 
and n are real constants and n > 0, are ruled out for any negative value of a. 
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Analyses of type la supernovae [H and of the cosmic 
microwave background, including the 3-year WMAP re- 
sults [2|, confirm that the present expansion of the uni- 
verse is accelerated. This is currently explained by in- 
voking a form of dark energy comprising 70% of the to- 
tal cosmic energy density p and with exotic properties (it 
does not cluster at small scales and has negative pressure 
Pde — —pde)- Recently, as an alternative to this exotic 
dark energy, it has been proposed to explain the cosmic 
acceleration with geometry by modifying Einstein gravity 
at the largest scales (or at low curvatures) by introducing 
corrections to the Einstein-Hilbert Lagrangian ( u f(R)", 
or "fourth order" , or "modified" gravity) . The simplest 
form for the action of modified gravity is 



S : 



-gf(R) 



(1) 



where K = 8ttG, R is the Ricci curvature [3j, and 
is the matter action. One can vary the action with re- 
spect to the metric ("metric formalism") or perform a 
Palatini variation in which the metric and the connec- 
tion are treated as independent variables Q ("Palatini 
formalism") — the resulting field equations, which coin- 
cide in general relativity, are different in higher order 
gravity. Furthermore, if the matter action S^ m ' also de- 
pends on the connection, one obtains a third possibility, 
metric-affine gravity theories . 

Quadratic quantum corrections were originally intro- 
duced to improve the renormalizability of general rela- 
tivity [7[ and were used to achieve inflation in the early 
universe Q . There is also motivation for modified gravity 
from string/M theory 

The field equations obtained by varying the action (Q]) 
in the metric formalism are 



f'Rab - ^ 9ab = V Q V 6 /' 



(2) 



where a prime denotes differentiation with respect to R, 
g a b and i? a & are the metric and Ricci tensors, and T a b 
is the stress-energy tensor of ordinary matter. The idea 
of modified gravity consists in introducing corrections to 
the Einstein-Hilbert action that are negligible in the early 



universe and only become effective at late times when the 
Ricci curvature R decreases (the radiation era with R = 
deserves a special discussion — see below). The proto- 
type is the theory f(R) = R- p A / R JrUS, El which, 
however, was found to be unstable 141 1 . Therefore, 
we parametrize the deviations from Einstein gravity as 

f(R)=R + e<p(R), (3) 

where e is a small parameter with the dimensions of an 
inverse squared length and (p is arranged to be dimen- 
sionless (in the previous example with f = R — p / R it 
is e = —p A with p ~ Ho w 10~ 33 eV, the present value 
of the Hubble parameter H). In order for the effects of 
these corrections to appear only at recent times a sen- 
sible choice is e ~ Hq ~ 10~ 66 eV 2 , which amounts to 
some fine-tuning which is present in all these theories, 
as well as in dark energy models. At present there is 
little indication on what the correct modification of Ein- 
stein gravity should be and it would be useful to weed 
out some of the competing models. Besides compatibil- 
ity with experimental data (e.g., [16l| ) . minimal criteria 
that a modified gravity theory must satisfy in order to be 
viable are the following: 1) reproducing the desired dy- 
namics of the universe including an inflationary era, fol- 
lowed by a radiation era (required and well-constrained 
by primordial baryogenesis and nucleosynthesis 17] and a 
matter era (required for the formation of structures) and, 
finally, by the present acceleration epoch (constrained 
by cosmic microwave background experiments, supernova 
data, and large scale structure surveys). There is re- 
cent evidence HH (see also [l9j]) that the popular models 



f(R) = R - p 2{ - n+1 )/R n with fi,n > |10|, lUj, l3J| do 
not allow a matter-dominated era and are therefore ruled 
out. 2) Having Newtonian and post-Newtonian limits 
compatible with the available Solar System experiments. 
This issue is not settled in general, there are conflict- 
ing results on the weak-field limit of modified gravity in 
specific scenarios, and we still lack a general formalism 
applicable to all modified gravity models 2(| 21 , 23] ■ 3) 
The theory must be stable at the classical and quantum 
level. There are in principle several kinds of instabilities 
to consider [l^ . Dolgov and Kawasaki [3] discovered 
an instability in the matter sector for the specific model 



2 



f(R) = i? — /i 4 /i?, a result confirmed in Ref. 14|, in 
which it is also shown that adding a term quadratic in 
R removes the instability. Instabilities of de Sitter s pac e 
in the gravity sector were found in Refs. 23j, [2J, |25|. l2q. 
while stability with respect to black hole nucleation was 
studied in Refs. 27j. The consideration of physically 



different instabilities yields remarkably similar stability 
conditions. The theory should also be ghost-free: the 
presence of ghosts has been studied in Refs. [2l|. In this 
communication we focus on the Dolgow-Kawasaki insta- 
bility and generalize their analysis to arbitrary f(R) the- 
ories. 

By taking the trace of eq. ([2|) one obtains 



3nf + f'R-2f = kT . 
By evaluating □/' this equation is written as 



tp'" (etp' - 1) 

aR+?—\7 a Rv a R+ v ; > R 

tp" Se tp" 



kT tp 
3etp" + 3p" 



(4) 



(5) 



We assume that the second derivative tp" is different from 
zero; if tp" — on an interval then the theory reduces to 
general relativity and is not of interest here. Isolated ze- 
ros of if" (at which the theory is "instantaneously general 
relativity") are not dealt with in this paper, but they are 
generally allowed and the following discussion applies to 
the distinct regions in which tp" maintains its sign. 

We now follow the procedure of Ref. [3] and consider 
a small region of spacetime in the weak-field regime, in 
which the metric and the curvature can locally be ap- 
proximated by 



9ab — Vab + h a0 , 



R = -kT + Rx 



(6) 



where r\ a b is the Minkowski metric and \R\/kT\ <C 1. 
This inequality excludes the case of conformally invariant 
matter with T = 0, including the important case of a 
radiation fluid with equation of state P = p/3 — this 
situation is examined later. Equation ([5]) yields, to first 
order in R\, 

9 2k tp" 1 ■ ■ 2k tp'" - 
i?i - V 2 i?i %— TRi + — %— VT ■ Vi?i 



Y 



Y 



(7) 



where V and V 2 denote the gradient and Laplacian oper- 
ators in Euclidean three-dimensional space, respectively, 
and an overdot denotes differentiation with respect to 
time. The function tp and its derivatives are now evalu- 
ated at R — — k T. The coefficient of R\ in the fifth term 
on the left hand side is the square of an effective mass and 
is dominated by the term (3e ip") 1 due to the extremely 
small value of e needed for these theories to reproduce 



the correct cosmological dynamics. It is therefore obvi- 
ous that the theory will be stable if ip" = f" > 0, while 
an instability arises if this effective mass is negative, i.e., 
if tp" = f" < 0. The time scale for this instability to 
manifest is estimated to be of order 10 -26 s for the spe- 
cific case <p{R) = —\x^jR 13]. As shown by this example, 
the smallness of <p" contributes to the magnitude of the 
effective mass and the smallness of the characteristic time 
scale for the instability to develop. Although a (p" term 
also appears in the denominator of the last source term 
on the right hand side of eq. , it is partially offset by 
the presence of tp and tp' in the numerator, while this 
does not occurs for the dominant term in the effective 
mass of Ri. 

Let us consider as an example the model f(R) = 
; — with n > and n > (n is not restricted to 



R- ^ 

be an integer): since /" = — n(n+ 1) ^ n+r > / R n+2 < 
this theory always suffers from the Dolgov-Kawasaki in- 



stability. Instead, the theory f(R) = R 



(again 



with /i, n > 0) is stable (but the required cosmological 
dynamics may not be achieved). As another example 
consider f(R) = R + A + aR n with n > and A > 
a cosmological constant, corresponding to Starobinsky 
inflation in the early universe [8j- The R 2 correction to 
the Einstein-Hilbert Lagrangian is large during this early 
epoch and the contribution of matter to the dynamics is 
negligible; any instability could in principle be interest- 
ing, and /" = naR n ~ 2 > if a > 0, corresponding to 
stability. 

We can now consider the case of a radiation fluid or 
other form of matter with vanishing trace T of the stress- 
energy tensor. In this case eq. ([7]) becomes 



Ri 



tp" 



i?i -V 2 i?i 



tp" 



V-Ri 



1 

"v 7 



- - <p 

e 



i?i = 



3tp" 



(8) 



Again, the effective mass term is w (2>etp") , which has 
the sign of /" and the previous stability criterion is re- 
covered [28jj. 

One could think of using an independent approach and 
study the equilibrium non-perturbatively by taking ad- 
vantage of the dynamical equivalence between f(R) grav- 
ity and a scalar-tensor theory. By introducing the auxil- 
iary field </>, the action ([1]) can be rewritten as 

' f d 4 xy/^ [iP(<j))R-V{(j))}+ S {m) (9) 



2k 

when /" ^ 0, where 

V#) = f'{4>) , 



v{<t>) = </>/'(</>) - m . (io) 



The corresponding field equations are 
1 / V 

Gab = — V a V b -0 - gabO^J - —g a b 

ip V 2 



K 

— Tab 

Ip 



(11) 
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R 



# dV 







(12) 



The action © reduces to ([T]) trivially if = R and, vice- 
versa, variation of @ with respect to 4> yields = i? if 
/" 7^ 0. Now one may think of studying the stability of 
the model by looking at the shape and extrema of the 
effective potential V(4>) but this would be misleading. In 
fact, the dynamics of <p are not regulated by V but are 
fixed by the strong constraint <fr = R. By taking the trace 
of eq. (fTTj) and using eq. (flUf . one obtains again eq. Q. 
Equation (fT2")l for (j> has no dynamical content because, 
using eq. (fT0|) with <^> = R, it is identically satisfied. This 
is linked to the fact that the kinetic term for <fi is absent 
in the action (J9j) with Brans-Dicke parameter oj — and 
the dynamical equation for (f> degenerates into an identity. 
Similar considerations apply to the Palatini action, which 
is equivalent to a Brans-Dicke theory with parameter u> = 
-3/2. 

In conjunction with the other theoretical requirements 
listed above, the (in)stability condition derived here con- 
stitutes a useful criterion to veto some of the proposed 
f(R) gravity scenarios from a theoretician's point of view, 
in order to focus on more promising models. 

Note that, when /" < 0, the instability of these theo- 
ries can be interpreted, following eq. (|7|), as an instability 
in the gravity sector. Equivalently, by using eq. ([6|), it 
can be seen as a matter instability (this is the interpreta- 
tion taken in [3]). Whether the instability arises in the 
gravity or matter sector seems to be a matter of inter- 
pretation. 

A problem potentially related to the instability is the 
weak-field limit of f(R) gravity, which has been the sub- 
ject of several papers reporting conflicting (sometimes 



even opposite) results [2CJ, |2l|, |22j. We hope to report 



soon on this issue and its relation with local instability. 

Finally, we mention the issue of the Cauchy problem 
a well posed Cauchy problem is necessary if the thcor 
is to have predictive power Q . It was shown in Ref. [2 
using harmonic coordinates, that theories of the form 



S = 



-g (R + aR ab R ah + f3R 



(13) 



have a well posed initial value problem. Moreover, the 
existence of a well posed Cauchy problem can be re- 
duced to the analogous problem by using the dynami- 
cal equivalence with scalar-tensor gravity © and (fTU|) 
when /" ^ 0. The well posedness of the Cauchy problem 
was demonstrated for particular scalar-tensor theories in 
Refs. [3(J[2l|] and has recently been the subject of a thor- 
ough analysis [23] . Although further study is needed for 
the specific cases u — and w = —3/2, which are rele- 
vant for f(R) gravity in the metric and Palatini formal- 
ism, respectively, it seems that the Cauchy problem is 



well posed [32] . Details on this subject will be presented 
elsewhere. 

The author thanks the organizers of the 2006 Mod- 
ern Cosmology Workshop in Benasque and the Benasque 
Center for Science, where this work was carried out, T.P 
Sotiriou for pointing out an error in a previous version 
of this paper, and financial support from the Natural 
Sciences and Engineering Research Council of Canada 
(NSERC). 



[1] 



[2] 
[3] 



[6] 
[7] 



[10] 

[11] 
[12] 

[13] 
[14] 



vfaraoni@ubishops.ca 

A.G. Riess et al., Astron. J. 116, 1009 (1998); Astron. 
J. 118, 2668 (1999); Astrophys. J. 560, 49 (2001); As- 
trophys. J. 607, 665 (2004); S. Perlmutter et al., Nature 
391, 51 (1998); Astrophys. J. 517, 565 (1999); J.L. Tonry 
et al, Astrophys. J. 594, 1 (2003); R. Knop et al, As- 
trophys. J. 598, 102 (2003); B. Barris et al, Astrophys. 
J. 602, 571 (2004). 



D.N. Spergel et al, |astro -ph/0603449 

In our notations, which follow Ref. 3, R = 6 ^ + 

in a spatially flat Friedmann-Robertson- Walker universe 

with scale factor a(t), an overdot denoting differentiation 

with respect to the comoving time t. 

R.M. Wald, General Relativity (Chicago University 

Press, Chicago, 1984). 

D.N. Vollick, Phys. Rev. D 68, 063510 (2003); Class. 
Quant. Grav. 21, 3813 (2004); X.H. Meng and P. Wang, 
Class. Quant. Grav. 20, 4949 (2003); Class. Quant. Grav. 
21, 951 (2004); Phys. Lett. 584B, 1 (2004); E.E. Flana- 
gan, Phys. Rev. Lett. 92, 071101 (2004); T. Koivisto, 
Phys. Rev. D 73, 083517 (2006); Class. Quant. Grav. 
23, 4289 (2006); T. Koivisto and H. Kurki-Suonio, Class. 
Quant. Grav. 23, 2355 (2006); P. Wang, G.M. Kremer, 
D.S.M. Alves, and X.H. Meng, Gen. Rel. Grav. 38, 517 
(2006); G. Allemandi, M. Capone, S. Capozziello, and 
M. Francaviglia, G en. Rel. Grav. 38 , 33 (2006); S. Nojiri 
and SD. Odintsov, fhep-th/0601213"l 
T.P. Sotiriou, Class. Qu ant. Grav. 23, 5117 (2006); T.P. 
Sotiriou and S. Liberati, gr-qc/0604006 ; N.J. Poplawski, 
Class. Quant. Grav. 23, 2011 (2006); 23, 4819 (2006). 
K.S. Stelle, Phys. Rev. D 16, 953 (1977); I.L. Buch- 
binder, S.D. Odintsov, and I.L. Shapiro, Effective Action 
in Quantum Gravity (IOP, Bristol, 1992). 
A. A. Starobinsky, Phys. Lett. 91B, 99 (1980). 
S. Nojiri and S.D. Odintsov, Phys. Lett. 576B, 5 (2003); 
S. Nojiri, S.D. Odintsov, and M. Sami, Phys. Rev. D 74, 
046004 (2006). 

S. Capozziello, S. Carloni, and A. Troisi, 
|astro-ph/0303041| 

S.M. Carroll, V. Duvvuri, M. Trodden, and M.S. Turner, 
Phys. Rev. D 70, 043528 (2004). 

X. Meng and P. Wang, Class. Quant. Grav. 20, 4949 
(2003); S.M. Carroll, A. De Felice, V. Duvvuri, D.A. 
Easson, M. Trodden, and M.S. Turner, Phys. Rev. D 71, 
063513 (2005). 

A.D. Dolgov and M. Kawasaki, Phys. Lett. 573B, 1 
(2003). 

S. Nojiri and S.D. Odintsov, Phys. Rev. D 68, 123512 
(2003); Gen. Rel. Gravity 36, 1765 (2004). 



4 



[15] J.D. Barrow and A.C. Ottewill, J. Phys. A 16, 2757 
(1983); A. Nunez and S. Solganik, Phys. Lett. 608B, 
189 (2005); |hep-th/0403T59} T. Chiba, J. Cosmol. As- 
tropart. Phys. 0503, 008 (2005); P. Wang, Phys. Rev. D 
72, 024030 (200 5); A. De Fe lice, M. Hindmarsh, and M. 
Trodden, |astro-ph/0604154| 

[16] O. Mena, J. Santiago, and J. Weller, Phys. Rev. Lett. 96, 
041103 (2006). 

[17] G. Lambiase and G. Scarpetta, Phys. Rev. D 74, 087504 
(2006). 

[18] L. Amendola, D. Polarski, and S. Tsujikawa, 
|astro-ph/0603703| 

[19] S. Capozziello, S. Nojiri, S.D. Odintsov, and A. Troisi, 
astro-ph/060443; L. Amendola, D. Polarski, and S. Tsu- 
jikawa, astro-ph/060443; S. Nojiri and S.D. Odintsov, 
|hep-th/0608008| A.W . Brookfield, C. van de Bruck, and 
L.M.H. Hall, |hep-th/0608015| 

[20] M.E. Soussa and R.P. Woodard, Gen. Rel. Grav. 36, 
855 (2004); R. Dick, Gen. Rel. Grav. 36, 217 (2004); 
A.E. Dominguez and D.E. Barraco, Phys. Rev. D 70, 
043505 (2004); D.A. Easson, Int. J. Mod. Phys. A 19, 
5343 (2004); G.J. Olmo, Phys. Rev. L ett. 95, 261102 

(2005) ; Phys. Rev. D 72, 083505 (2005); [gr^qc/0505135| 
gr-qc/0505136 ; I. Navarro and K. Van Acoleyen, Phys. 
Lett. 622B, 1 (2005); G. Allemandi, M. Francaviglia, 
M.L. Ruggiero, and A. Tartaglia, Gen. Rel. Grav. 37, 
1891 (2005); J.A.R. Cembranos, Phys. Rev. D 73, 064029 

(2006) ; S. Capozziello and A. Troisi, Phys. Rev. D 72, 
044022 (2005); T. Clifton and J.D. Barrow, Phys. Rev. 
D 72, 103005 (2005); T.P Sotiriou, Gen. Rel. Grav. 38, 
1407 (2006); C.-G. Shao, R.-G. Cai, B. Wang, and R.- 
K. Su, Phys. Lett. B 633, 164 (2006); S. Capozziello, A. 



Stabile and A. Troisi, |gr-qc/060307l| 
[21] I. Navarro and K. Acoleyen, J. Cosmol. Astropart. Phys. 

0603, 008 (2006); |gr-qc/0512109| G. Calcagni, B. de 

Carlos, and A. De Felice, Nucl. Phys. B 752, 404 

(2006); A. De Fel ice, M. Hindmarsh, and M. Trodden, 

|astro-ph/0604154l 
[22] V. Faraoni, Phys. Rev. D 74, 023529 (20 06); A. Ericke ck, 

T.L. Smith, and M. Kamionkowski, astro-ph/0610661 , G. 

Allemandi and M.L. Ruggiero, astro-ph/0610661 
[23] V. Faraoni, Phys. Rev. D 70, 044037 (2004); Phys. Rev. 

D 72, 061501(R) (2005); V. Faraoni and S. Nadeau, Phys. 

Rev. D 72, 124005 (2005). 
[24] A. Dolgov and D.N. Pelliccia, Nucl. phys. B 734, 208 

(2006). 

[25] J.D. Barrow and S. Hervik, Phys. Rev. D 73, 023007 
(2006). 

[26] V. Miiller, H.-J. Schmidt, and A. A. Starobinsky, Phys. 

Lett. 202B, 198 (1988); H.-J. Schmidt, Class. Quant. 

Grav. 5, 233 (1988); A. Battaglia Mayer and H.-J. 

Schmidt, Class. Quant. Grav. 10, 2441 (1993). 
[27] G. Cognola, E. Elizalde, S. Nojiri, S.D. Odintsov, and S. 

Zerbini, J. Cosmol. Astropart. Phys. 02, 010 (2005); J. 

Phys. A 39, 6245 (2006). 
[28] A stability condition similar to /" > was obtained in 

Ref. |29j] in the metric formalism, and for cosmological 

perturbations in the Newtonian gauge. 
[29] Y.-S. Song, W. Hu, and I. Sawicky, [astro-ph/0610532" 
[30] D.R. Noakes, J. Math. Phys. 24, 1846 (1983). 
[31] W.J. Cocke and J.M. Cohen, J. Math. Phys. 9, 971 

(1968). 

[32] M. Salgado, Class. Quant. Grav. 23, 4719 (2006). 



